Vortices in a cylinder: Localization after depinning 
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Edge effects in the depinned phase of flux lines in hollow 
superconducting cylinder with columnar defects and electric 
current along the cylinder are investigated. Far from the ends 
of the cylinder vortices are distributed almost uniformly (de- 
localized). Nevertheless, near the edges these free vortices 
come closer together and form well resolved dense bunches. 
A semiclassical picture of this localization after depinning is 
described. For a large number of vortices their density p(x) 
has square root singularity at the border of the bunch (p(x) 
is semicircle in the simplest case). However, by tuning the 
strength of current, the various singular regimes for p(x) may 
be reached. Remarkably, this singular behaviour reproduces 
the phase transitions discussed during the past decade within 
the random matrix regularization of 2d-Gravity. 

PACS numbers: 74.60.Ge, 72.15.Rn, 05.45. +b 



Statistics of ensembles of flux lines in high-T c super- 
conductors has been a subject of numerous experimental 
and theoretical investigations (see for review Ref. [|). An 
important modification of the superconductor is achieved 
by creation of artificial disorder in the form of columnar 
pins produced by heavy ion irradiation. The vortex dis- 
tribution in the cylindrical superconductors with colum- 
nar defects and longitudinal electric current has attracted 
a renewed interest recently The current creates a 
transverse magnetic field, which attempts to wind vor- 
tices around the cylinder. For low current the vortices 
are trapped by the fluctuations in the density of defects 
and do not curl (transverse Meissner effect j| ) . With the 
increase of current the transition to depinned phase takes 
place. The mapping of this transition onto the derealiza- 
tion transition in Id non-Hermitean quantum mechanics 
Igl has caused an immediate and wide interest Q. How- 
ever, as it was pointed out in Ref. even after the tran- 
sition to complex spectrum eigenfunctions of correspond- 
ing non-Hcrmitcan Hamiltonian still exhibit the features 
of both localized and delocalized states. Physical conse- 
quences of this " localization after depinning" for vortices 
in the thin-shelled cylinder will be the subject of this 
paper. Our main prediction is illustrated by the Fig. 1. 
The depinned fraction of vortices is practically uniformly 
distributed over the surface far from the ends of cylin- 
der. However, while approaching the ends, vortices come 
closer and form well resolved localized bunches. The den- 
sity of vortices in this bunch at the edge coincides with 
the density of eigenvalues for the ensembles (Gaussian 
or non-Gaussian) of orthogonal random matrices. The 
phase transitions emerging while tuning the parameters 



in these random matrix models were a subject of inten- 
sive investigation in the past decade in the context of 
problems of 2d-Gravity H. As we show in this paper, 
experimental investigation of vortices in cylindrical sam- 
ples may open the way for the direct observation of such 
transitions in real low-dimensional systems. 

The classical energy for a flux line in the cylinder may 
be written as an action of equivalent particle H 
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The length and period of the cylinder are L T and I, 
x+l = x. The potential V(x) accounts for the interaction 
of vortex with the columnar defects and h is proportional 
to the longitudinal current. All our results are valid for 
random V(x). However, in order to enhance the effect it 
is better to prepare the sample with smooth and inho- 
mogeneous at the scale Ax ~ I density of columnar pins. 
The partition function for classical vortex now takes the 
form of quantum evolution operator in imaginary time 
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where [3 = 1/T is the inverse temperature and we have 
introduced the non-Hcrmitcan Hamiltonian 



H 



]_d_ 



V(x) 



(3) 



Thus the electric current in the cylinder acts like an imag- 
inary vector potential ih. The path integral in Eq. (|^) 
includes strings with given end points x\ and xi- In or- 
der to find the partition function for vortex with free ends 
one has to integrate over xi,2 0- The evolution operator 
Eq. (0) may be written in the form 



(4) 



where the left- and right- eigenvectors are defined via 
Hip R = eip R and H T ip L = eip L . The solutions are nor- 
malized via J Q ipiip^dx = 5ij. In this paper we will con- 
sider only the case h^ 1 <C I <C L T . In particular this 
means that only the ground state contribution survives 
in the sum in Eq. (Q). 

Interacting vortices in this approach are equivalent 
to the interacting bosons. Moreover, we will consider 
only the case of strong (^-function repulsion of vor- 
tices (bosons). 
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Interesting for us physical information is contained in 
the slow harmonics of the potential V(x). Therefore, 
in this paper we will restrict our attention on the case 
dV/dx <§; V . Also we will consider the case of strong 
depinning current h 2 3> V. 

exp(— a) 



a = / (y/2(V - e) - h)dx' + ... . (5) 
Jo 



The corrections to a may be also easily found. The quan- 
tization condition, which allows one to find the set of 
complex eigenvalues e n , is a (I) — c(0) = 2nni. It will be 
enough for us to consider only low lying excitations of 
the Hamiltonian Eq. (g). In this case 

e n = -(ikT - hf/2 + {V) , k = 2nn/l , (6) 

r* V(x')-(V) 
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Here (V) = l^ 1 J Q Vdx. The most exciting in this re- 
sult is that <j{x) acquires the nontrivial real part in the 
presence of columnar field V(x) (in general \Rea\ ^ 1 if 
I ^S> h/V). This means that even after the depinning 
transition eigenvectors of the non-Hcrmitcan Hamilto- 
nian Eq. (^) still remain strongly (exponentially!) lo- 
calized. The localization described by the Eqs. ((5||) 
has almost nothing in common with usual localization 
at h = 0. To see this consider the simple example of 
very wide square- well potential: V = Uq for < x < 1/2 
and V = —Uo for — 1/2 < x < 0. Instead of almost 
extended ground state in the absence of transverse field 
ijj(h = 0) ~ sin(27rir/Z) at — 1/2 < x < one finds via the 
Eq. (p) the clearly localized even for h 2 ^> Uq eigenvec- 
tor tpQ ~ exp(— Uo\x\/h) ||. Moreover, we see from this 
example that the imaginary vector potential first creates 
the localized states (maximum of localization is reached 
at h 2 ~ V — (V)). However, with further increase of h 
the localization length grows up again. All eigenstates 
described by the Eq. (^) are localized near the minimum 
of ao(x). This feature of " non-Hermitean" localization 
is also in sharp contrast with the usual Anderson case, 
where eigenstates with close energy strongly repel in the 
coordinate space. 

The Eqs. (pl[|) describe localization of eigenvectors 
by the long wave-length harmonics of potential V(x) in 
strong imaginary vector potential. In the case of ran- 
dom disorder, however, a part of eigenfunctions (whose 
with the localization length £ < h^ 1 ) may be localized in 
the usual Anderson sense on the rare local fluctuations 
of V(x) M. In this case, if there are only few vortices in 
the sample, they will be trapped by these true localized 
states. However, with increase of the number of fluxes, 
the fraction of depinned vortices emerges. Correspond- 
ing wave functions are described again by the analog of 
Eqs. (||[3|). Just this situation is shown on the Fig. 1. 

The localization described by the WKB formula 
Eq. (|s|) should have a simple classical explanation. In- 
deed, the action for the solution of classical equation of 



motion for the string connecting points X\ and X2 with n 
windings over the cylinder may be shown to have a form 
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where the classical energy e should be found from 
Ix! +X2 dx/ y/2(V — e) = L T . Note that we use the re- 
versed sign of the kinetic energy e = — i 2 /2 + V com- 
pared to usual classical definition. In order to find the 
most important classical configuration we have to deter- 
mine the minimum of S^ n ' with respect to n. For n» 1 



one has for A5 C 
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AS cl = Se (L T -L T )+ (V2(F - e) - h)dx , (8) 



where Se is the energy difference for two paths. Thus, the 
condition of extremum dS^ /dn = coincides with the 
semi-classical quantization condition for the ground state 
of Hamiltonian Eq. (^). Consequently, the contribution 



(n) 

of smallest into the functional integral _Eq 



is 



enough to reproduce the WKB result Eqs. (D,|[) in the 
limit of large L T . The classical consideration allows also 
to find the typical tilt of the vortex trajectory (x) — 
h + 0(V 2 /h 3 ). 

The Hamiltonian Eq. (|j|) and its transpose are re- 
lated via H T (h) = H{-h) @. This simple equal- 
ity allows one to introduce the conserving current for 
left- and right- eigenvectors with the same energy J = 
iP L V fl - ip L ip R ' + 2hiP L ^ R ; J' ee dJ/dx = 0. For small 
h and random V(x) in the thermodynamic limit I — * oo 
one has J = for all low (real)energy states. Above 
some critical value Res = e c the current become non- 
zero and the energy acquires nontrivial imaginary part. 
This appearance of conserving current J ^ was usu- 
ally considered as an indication of existence of delocal- 
ization transition in Id non-Hermitean quantum mechan- 
ics. However, formally the existence of current shows the 
absence of exponential localization only for quantities bi- 
linear in ip L and ip R - Individually ip L and ip R still may 
be localized. Indeed, for the WKB wave functions Eq. (|^) 
one finds ip L (x)ip R (x) — l" 1 , i.e. exponential growth of 
i/j r is compensated by the same decrease of ip L . The 
probability to find vortex at the point x on the trans- 
verse slice r of the cylinder due to the Eq. (0) has the 
former) = Z- 1 {^f\e- p{ - L ^ T ^ H \x){x\e-' fiT "\^ 1 ). For 
large r, L T — t ~ L T one finds p(x,r) ~ ip R (x)tpQ (x). 
We see that the effect of localization is washed out in- 
side the cylinder. However at the top of cylin- 
der p R (x) — p(x,L r ) ~ iPq{x) and at the bottom 
p L {x) = p(x,0) ~ tpo(x) the vortex is strongly localized. 

As we have told, the ensemble of vortices should be 
treated as the system of interacting bosons with single 
particle Hamiltonian Eq. (|^). At least for low density of 
fluxes, the inter-boson interaction is equivalent to short- 
range repulsion U — UqS(x — x'). For many vortices the 
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p r 

single-particle wave functions ipi ' and energy £j should 
be replaced by the many-particle ones ^ i ' and the total 
energy E. The ground state for impenetrable bosons is 
found by the analogy with iV-fermion solution [nof 
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Here ip R are the eigenvectors of the single particle Hamil- 
tonian Eq. (||) with lowest Ree n . For odd N the eigen- 
vectors should be the usual periodic ip R {x + = *l J n( x )- 
However, for even TV one should use the antiperiodic 
boundary conditions ipn( x + = ~i } ni x )- ^ n both cases 
the total ground state energy E = ^2 e n entering the 
analog of the Eq. (Q) is real. Also for the even and odd 
N and smooth V(x) the 4>n( x ) are we -^ described by the 
Eqs. (^|J^)- Within this approximation the wave function 
is further simplified 



n 



11* 



(10) 



Furthermore, in the most interesting case, when the func- 
tion exp(— cto(x to )) has a narrow (Ax -C I) peak one may 
replace | sin(7r(xi — Xj)/l)\ by ir\xi — Xj\/l. This means 
that the distribution of vortices described by the ^> R co- 
incides formally with the distribution of eigenvalues of 
N x N Orthogonal Random Matrix Ensemble with the 
weight function exp(— Trao(M)). The density of vortices 
p(x) at the end of the cylinder is found after integration 
of the ^ R over the all Xi except for one. The saddle-point 
method (large N approximation) for calculation of such 
integrals was developed many years ago in Ref. pd[ . In 
particular, in the most general case of quadratic mini- 
mum of <tq one has 



a(x — xq) 2 
a ° = 2hT ; 9 
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Here we have shown explicitly the dependence of <jq on 
the transverse field h and temperature T. The coefficient 
a is determined only by the potential V(x). With the 
increase of N the anharmonic contributions to a should 
be also taken into account. For example one may write 
(we put xq = 0) 



1 ax 2 



W{x/X) , W(0) = 1 . 



(12) 



Now a contains the information about the strength of in- 
teraction V(x), while A is the characteristic length for its 
variation. The vortices at the edge of the cylinder in this 
case also form a smooth dense bunch with square root 
p(x) at the border (see Fig. 2a). The width of the bunch 
is x c ~ \jNKT/a and the anharmonic contributions be- 
came important starting from N ~ N c — a\ 2 /hT. The 
new phenomena may take place if the function ctq(x) has 



more than one minimum. Some variants of a peculiar 
behavior of p(x) in this case are illustrated by Figs. 2b 
and 2c. With increase of N in the case of two minimums 
the second small stable bunch of vortices is born at some 
value N\. These two bunches join together at the second 
critical value N2- The moment of consolidation of two 
bunches into one is shown on the Fig. 2b. The density 
of states close to the transition in the vicinity of meeting 
point xq is 



p = a(x — xq)\/ (x — xo) 2 + A[h — h c ] for h < h c , 

p = a((x - x a ) 2 + A[h - h c }/2) for h> h c , (13) 

where a, A, xq and h c vary smoothly with the change 
of N, T, or h. The number of vortices which may be 
kept in equilibrium in each well may be regulated by the 
external parameter h. This feature of the vortex distri- 
bution opens the way to create the metastable configu- 
rations like one shown on the Fig. 2c. For h < h* the 
vortices are confined (during exponentially long time) in 
the deepest well. Above h* the fast decay into second 
well takes place. The critical configuration shown on the 
Fig. 2c is characterized by the novel singular behaviour 
of the density at the border p(x) ~ (x — xq) 3 ^ 2 . Just this 
type of critical behaviour corresponds to the continuous 
limit in the Random Matrix regularization of 2d-Gravity. 

The catastrophic change in vortex distribution at the 
transitions should also change the thermodynamic char- 
acteristics of the ensemble of flux lines. After integration 
over the positions of the ends of vortices one finds the 
partition function and the Free energy 



F = -Tln(Z) = F Lt 



F R + Ft 



(14) 



The contribution proportional to the total volume Fl t 
in our simple case has the form 

F Lt =NL t {-h 2 /2 +(V) + (ttNT/1) 2 /6 + /„} . (15) 

Here fo(T) accounts for the short wave-length fluctua- 
tions of string and does not depend on h, V(x) and N. 
The more interesting for us are the edge contributions 
F R and F^ . For large N and the quadratic ctq 



3 1, / a R l 2 \ , 



Here again go depends on the details of regularization of 
the functional integral Eq. (0). With the increase of N 
(or h, T) the edge Free energy changes smoothly until one 
meets one of the singular points considered above. The 
first possible singularity is the birth of new small bunch. 
The corresponding correction at N > Ni is 



AF e = -AT(N-Nif 



(17) 



Here A = A(N, hT) > 0. As we have learned from the 
2d-Gravity, this correction is purely nonperturbative, i.e. 
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it could not be related with the analytic behaviour of 
F e below the singularity. Another two kinds of singular 
corrections associated with the confluence of two bunches 
(h c ) and decay of metastable bunch (h*) lead to 



AF e 
N 2 T 



(h c -h) 3/2 and 



AF* 
N 2 T 



(K - h) 



5/2 



(18) 



The constants h c and h* are the functions of N and T. 
Depending on the concrete way of realization of physical 
experiment one may write instead of the Eq. ([fq), for 
example, AF - (T c - T) 3 / 2 or AF ~ (N 2 - Nfi*. 

In summary, the edge effects in vortex distribution in 
superconducting cylinders may provide us with the vari- 
ety of new phenomena with clear experimental signature 
(see again the Figs. 1,2). Among them are the strong 
localization of vortices at the end of cylinder and vari- 
ous critical regimes for this localization available by tun- 
ing of the longitudinal current. Technically these effects 
arose due to the peculiar features of localization in non- 
Hermitean quantum mechanical Hamiltonian d§|). From 
the pure theoretical point of view, the most exciting is 
the correspondence between distribution of flux lines at 
the end of cylinder and distribution of eigenvalues of the 
ensembles of random matrices. The ensemble of fluxes 
turns out to be the almost unique example of the system, 
where not only local (correlations of close levels etc.), but 
also global features of random-matrix spectrum are of 
100% importance. For example, the phase transitions in 
non-Gaussian Matrix ensembles have been a subject of 
enormous activity in last 10 years within the context of 
2ci-Gravity ||. However, to the best of my knowledge, in 
this paper the first proposal is presented of a real phys- 
ical experiment where such a singular behaviour may be 
observed. 
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FIG. 1. The vortex distribution on the cylinder (I = 0). 
The Hamiltonian H has one eigenstate localized in the usual 
sense due to a strong local fluctuation of V(x). In the typical 
transverse section far from the ends of the cylinder one vortex 
is trapped by this localized level, while others are almost uni- 
formly distributed. However, near the ends these free vortices 
come closer together and form well resolved dense bunches. 
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FIG. 2. Examples of distribution of vortices at the end of 
cylinder: a). Semicircle in a single well b). Confluence of two 
bunches in a double well c). Opening of a decay of metastable 
bunch. 
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